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Abstract 



We have carried out extensive series studies, at T = and at high tempera- 
tures, of 2-chain and 3-chain spin-half ladder systems with antiferromagnetic 
intrachain and both antiferromagnetic and ferromagnetic interchain couplings. 
Our results confirm the existence of a gap in the 2-chain Heisenberg ladders 
for ah non-zero values of the interchain couphngs. Complete dispersion rela- 
tions for the spin-wave excitations are computed. For 3-chain systems, our 
results are consistent with a gapless spectrum. We also calculate the uniform 
magnetic susceptibility and specific heat as a function of temperature. We 
find that as T ^ 0, for the 2-chain system the uniform susceptibility goes 
rapidly to zero, whereas for the 3-chain system it approaches a finite value. 
These results are compared in detail with previous studies of finite systems. 

PACS Indices: 75.10.-b., 75.10J., 75.40.Gb 
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I. INTRODUCTION 



The magnetic properties of low dimensional systems have been the subject of intense 
theoretical and experimental research in recent years. It is by now well established that 
one-dimensional Heisenberg antiferromagnets with integer spin have a gap in the excitation 
spectrum, whereas those with half-integer spin have gapless excitations. The former have 
a finite correlation length, while for the latter it is infinite with the spin-spin correlation 
function decaying to zero as a power law. In 2-dimensions, the unfrustrated square-lattice 
Heisenberg model has long range Neel order in the ground state. It has gapless Goldstone 
modes as expected. In recent years much interest has focussed on systems with intermediate 
dimensionality and on questions of crossovers between d = 1 and d = 2. One approach 
to this problem has been to study a two-dimensional system where the coupling for the 
spins separated along the x-axis is different from those separated along the y-axis flj. It 
has been suggested that an alternative way to explore this issue is through the Heisenberg 
spin ladders consisting of a finite number of chains coupled together, with a coupling J^^ 
along the chains and Jj_ between them. These systems have been the subject of considerable 
recent theoretical and experimental interest. 

Experimentally, 2-chain 5" = | ladders are realized in vanadyl pyrophosphate (^0)2^207 
and in the strontium cuprate SrCu20^ 0, whereas 3-chain S = \ ladders are realized 
in the strontium cuprate SrCu^O^ 0. 

Theoretically, a number of striking predictions have been made for such systems. These 
have been recently reviewed by Dagotto and Rice . Barnes et al. |^,^ were the first to carry 
out extensive Monte Carlo studies of the excitation spectrum and the magnetic susceptibility 
for 2-chain ladders with antiferromagnetic interchain coupling. White [0] et al. and Hida 
have used the density matrix renormalization method to study the spin gap. Watanabe |^ 
has applied the numerical diagonalization method to finite systems of 2-chain ladders with 
ferromagnetic interchain coupling. Azzouz et al. developed a mean-field theory and used 
the density-matrix renormalization group method to study 2-chain ladders. Gopalan, Rice 
and Sirgrist [11| presented a variational wavefunction for the ground state of the 2- leg ladder. 
Rojo [0 has shown that spin ladders with odd number of legs have gapless excitations. 
Troyer et al. [|13[ have used improved versions of the quantum transfer-matrix algorithm to 
study the temperature dependence of the susceptibility, specific heat, correlation length etc 
of 2-chain ladders. Finite-size scaling was used by Hatano for multi-leg ladders, and 



recently Frischmuth et al. |]15[ and Sandvik et al. ||I6[ have applied Quantum Monte Carlo 
simulation to compute the temperature dependence of uniform susceptibility and internal 
energy for spin ladders with up to 6 legs. One clear result from all these studies is that 
ladders with an even number of legs have an energy gap, short range correlations and a 
"spin liquid" ground state. On the other hand, ladders with odd number of legs have 
gapless excitations, quasi long range order, and a power-law falloff of spin-spin correlations, 
similar to single chains. Experiments also confirm these features. 

We have carried out extensive series studies of 2-chain and 3-chain ladder systems with 
both antiferromagnetic and ferromagnetic interchain coupling J±, via Ising expansions and 
dimer expansions at T = 0, and also by high temperature series expansions. Our results 
confirm the existence of a gap in the 2-chain system and delineate the phase diagram in the 
parameter space of Ising anisotropy and the parameter ratio J±/J,,. The complete spin- 
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wave excitation spectra are computed. For the 3-chain system we are unable to exclude the 
possibility of a small gap, although our results are consistent with a gapless spectrum. In 
addition, we develop a high-temperature series expansion for the uniform magnetic suscepti- 
bility and the specific heat for 2-chain and 3-chain systems with J = J±; the susceptibility 
of 2-chain ladders is as expected for a system with a spin-gap while that of 3-chain ladders 
appear to remain finite in the zero-temperature limit suggesting an absence of a spin-gap. 
We compare our results in detail with previous calculations. 



II. SERIES EXPANSIONS 

The Hamiltonian of a spin ladder with rii legs is given by, 

l=ni l=ni — l 

H = J // Si^i ■ S;,i+i + J_L ^ Si^i ■ Sj+i^i (1) 
i,l=l i,l=l 

where S;^j denotes the S=l/2 spin at the ith site of the /th chain. J^^ is the interaction 
between nearest neighbor spins along the chain and Jj_ is the interactions between nearest 
neighbor spins along the rungs. We denote the ratio of couplings as y, that is, y = J_l/ J^^. In 
the present paper the intrachain coupling is taken to be antiferromagnetic (that is, J^^ > 0) 
whereas the interchain coupling J± can be either antiferromagnetic or ferromagnetic. This 
includes the values of interest in the experimental systems discussed earlier where J± ~ J^^ . 
Without loss of generality, we can set J^^ = 1 hereafter. 

We have carried out three different expansions for the system. The first is the expansion 
about the Ising limit at zero temperature for both two and three-chain ladders. We have 
computed the ground state properties as well as the spin-wave excitation spectra by this 
expansion. The second is the dimer expansion, again at T = 0. This expansion can be done 
for the 2-chain system with antiferromagnetic interchain coupling only. The third is the 
high-temperature series expansion for the uniform susceptibility of the 2-chain and 3-chain 
ladders with y = 1. 



A. Ising expansions 

To perform an expansion about the Ising limit for this system, we introduce an anisotropy 
parameter x, and write the Hamiltonian in Eq.(|l|) as: 



where 



H = Ho + xV (2) 



l=ni l=ni — l 

i,l=l i,l=l 

l=ni l=ni — l 

v=y: + ^M^M+i) + y E (stiSf+i, + Si.sUi.) (3) 

i,l=l i,l=l 
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The limits x = and x = 1 correspond to the Ising model and the Heisenberg model re- 
spectively. The operator Hq is taken as the unperturbed Hamiltonian, with the unperturbed 
ground state being the usual Neel state for antiferromagnetic interchain coupling, and a fully 
ordered collinear state for ferromagnetic interchain coupling. The operator V is treated as 
a perturbation. It flips a pair of spins on neighbouring sites. The Ising expansion method 



has been previously reviewed in several articles [|17|p!8| , and will not be repeated here. The 
calculations involved a list of 9184 linked clusters of up to 16 sites for the 2-chain ladder, 
and 14082 linked clusters of up to 12 sites for the 3-chain ladder, together with their lattice 
constants and embedding constants. 

Series have been calculated for the ground state energy per site Eq/N, the staggered 
magnetization M for y > (or collinear magnetization M for y < 0), the parallel stag- 
gered/colinear susceptibility X//^ the uniform perpendicular susceptibility x P^r site for 
several ratio of couplings y = ±0.1, ±0.25, ±0.5, ±0.75, ±1, ±1.5, ±2, ±4, ±8 up to order x^^ 
for 2-chain ladders, and for 3-chain ladders (the series for uniform perpendicular suscep- 
tibility X is one order less in each case). The resulting series for y = ±1 for the 2-chain and 
3-chain systems are listed in Tables I and II, the series for other value of y are available on 
request. 

To analyze these series, we first performed a standard Dlog Pade analysis of the mag- 
netization M and parallel susceptibility X// series. We found that for 2-chain ladders, the 
series lead to a simple power- law singularity at x < 1: 

M~(l-x/x,f x,-{l~x/x,)-^ (4) 

with the indices (3 and 7 close to 1/8 and 7/4 respectively. This transition at a; < 1, with 
criticality in the universality class of the 2Z)-Ising model, is strong evidence that in the 
Heisenberg limit the system has a disordered ground state with a spin-gap, as is the case 



for the spin-one chain |]I9[. In contrast for the 3-chain ladders, the series analysis showed 
poor convergence and suggested a singularity at x > 1. This implies that for the 3-chain 
ladders, the system is analogous to the spin-half chain, with gapless spectra and power-law 
correlations in the Heisenberg limit. 

Fig. 1 shows the phase boundary for 2-chain ladders as function of y. It is interesting 
to note the different behaviour for y > and for y < 0: for antiferromagnetically coupled 
2-chain ladders {y > 0), Xc decreases as y increases, and in the limit of y 00, Xc will 
approach 0. But for ferromagnetically coupled 2-chain ladders {y < 0), Xc first decreases as 
the absolute value of y increase from zero, but then the trend reverses and it, once again, 
approaches 1 as ?/ — — cxo. To understand this behavior, we can map the system for large 
negative y to a spin-one chain with on-site ( single-ion) anisotropy: 

H = l T.[^^,tot ■ S.+i,to* ± y{l - x){Sl,,,)'] (5) 
^ i 

where Si^tot denotes the S=l spin at the ith site of the chain. We can get the asymptotic 
behaviour of the phase boundary in the limit of y —00 hj studying the following spin-one 
chain with on-site anisotropy: 
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For this model, we have carried out series expansions in x' to order a;'^^ (i.e. 14 sites) 
for staggered magnetization M for several different value of c: c = 0.275, 0.29, 0.3, 0.325, 
and performed a standard Dlog Fade analysis to find the critical value c' which gives the 
singularity of M at x' = 1. We get d = 0.29(1). Hence, the asymptotic behaviour of the 
phase boundary in the limit ?/ ^ — oo is given by 

Xc=l + 0.29/y (7) 

which is also shown in Fig. 1 as a bold line near y/{l + \y\) = —1. 

Fig. 2 gives the results of the ground-state energy per site Eq/N versus y for both 2-chain 
and 3-chain ladders at the Heisenberg point x = 1. Our results for the ground state energy 



agree extremely well with the recent quantum Monte Carlo simulation [|1^]. In Fig. 3, we 
present the uniform perpendicular susceptibility at T = 0. 

We also performed the Ising expansion for the triplet spin-wave excitation spectrum of 
2-chain and 3-chain ladders using Gelfand's method pO|- To overcome a possible singularity 



at a; < 1 in the 2-chain ladders, and to get a better convergent series in the Heisenberg limit, 
we add the following staggered field term to the Hamiltonian in Eq. (0): 

^H = t{l-x)Y.{-irSt (8) 

i 

AH vanishes at x = 1. We adjust the coefficient t to get the most smooth terms in the 
series, with a typical value being t = 2. We computed the Ising expansion for the triplet 
spin-wave excitation spectrum e{k) up to order x^^ for 2-chain ladders, and up to order x^^ 
for 3-chain ladders. These series are too long to be listed here, but are available on request. 

These series have been analyzed by using integrated first-order inhomogeneous differential 
approximants For the 2-chain ladder, there are 2 bands of excitations. Fig. 4 shows 
the dispersion e{k), with ky = vr, for antiferromagnetic interchain coupling. The other band 
with ky = is related to this by e{kx, 0) = e(7r — k^, n). This is simply due to the staggered 
field, which doubles the spectrum. As a comparison, the dispersion relation of a single chain 
(that is the case of y = 0) is also shown. It can be seen from the graph that in the limit 
y ^ 0, the dispersion relation has a simple cosine function with a period of 27r, and in the 
limit y oo, the dispersion relation also has a simple cosine form with a period of Air, and 
a gap in the spectrum. Barnes and Riera argued that the dispersion, for all y, can be 
fitted by combining these two functions into the following form: 

e{k„ nf = e(0, rrf cos\kj2) + e(7r, uf sm\kj2) + cq sm\k,) (9) 

For y = 1, the Ising expansions give an energy gap of e(7r,7r) = 0.44(7). A more precise 
estimate is obtained by the dimer expansions, which give eiji, tt) = 0.504(7). We will discuss 
the dimer expansions later. 

For ferromagnetic interchain coupling, the two bands of spectra are independent, but 
each band is a simple cosine function with a gap at the minimum and symmetric about 
kx = 7r/2, as shown in Figs. 5-6. As noted in Figs. 4, it is clear that the spin-gap decreases 
smoothly as y is reduced, and vanishes at |/ = 0. These results agree well with previous 
calculations 0. 

For the 3-chain system, there are three bands. In the Ising limit, 2 bands have initial 
excitations located in the side rows, and the third band has it in the middle row. Figs. 7-12 
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show the spectrum of the three bands for ferromagnetic and antiferromagnetic interchain 
couphngs. From these graphs, we can see that all of the dispersion relations have a simple 
cosine function (except for middle row band with large y) with a minumum located at fc^; = 
(or fcj. = TT by symmetry), where two of these three bands have a definite gap, the third (the 
ky = band) is consistent with a gapless spectra. The estimate for the gap in the third 
band for all y values is 0.2(3) (except for the case of y = where we got 0.08(10)). Hence, 
we cannot exclude the possibility of a small gap here. But our results are consistent with a 
gapless spectrum, and given our earlier results on the phase boundary with Ising anisotropy, 
we believe the spectra are gapless. 



B. Dimer expansions 

For 2-chain ladders, with antiferromagnetic coupling between the chains, there is an 
alternative T = expansion that can be developed. In the limit that the exchange coupling 
along the rungs J_l is much larger than the coupling J^^ along the chains, that is ?/ 3> 1, 
the rungs interact only weakly with each other, and the dominant configuration in the 
ground state is the product state with the spin on each rung forming a spin singlet, so the 
Hamiltonian in Eq. (|I|) can be rewritten as, 

H/J^ = Ho + {l/y)V (10) 

where 

l=ni-l 
i,l=l 

(11) 

l=ni 
i,l=l 

We can treat the operator Hq as the unperturbed Hamiltonian. The eigenstates of a single 
pair of spins, or dimers, consists of one singlet state with total 5 = and eigen energy 
Es = -3/4: 

l^)s = ^(ITi)-UT)) (12) 

and three triplet states with total S = 1 and eigenenergy Et = 1/4: 

l*)* = [^(in) + UT)), ITT), lii)] (13) 

The operator V is treated as a perturbation. It can cause excitations on a pair of neigh- 
bouring dimers. Details of the dimer expansions and the matrix elements of V are given in 



Ref. [|18|, and will not be repeated here. 

We have carried out the dimer expansion for the ground state energy to order (l/y)^ 
and for the lowest lying triplet excitations to order (l/y)^. The series for the ground state 
energy per site Eq/N is: 
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Eo/N = J±[-3/8 - 3/(16y2) - 3/(32?/^) + 3/(256?/^) + 45/(512?/^) + 159/(2048?/^) 

-879/(32768y^) - 4527/(32768y^) - 248391/(2097152|/9)] (14) 

and the series for the excitation spectrum are hsted in Table III. Again, we use the integrated 
first-order inhomogeneous differential approximants to extrapolate the series. For y > 1, 



the dimer expansions give much better results than the Ising expansions. For y ~ 1 also 
the dimer expansions appear to converge better. The overall spectra determined from the 
combined study of dimer and Ising expansions are shown in Fig. 4. 



C. High temperature series expansions 

We now turn to the thermodynamic properties of the ladder system at finite temper- 
atures. We have developed high-temperature series expansions for the uniform magnetic 
susceptibility x{T) and the specific heat C{T), for 2-chain and 3-chain system with J± = J^^, 



C{T) 



df 



(15) 



where N is the number of sites, and /? = 1/(A;bT), and the internal energy U is defined by 

"-^^ (16) 

The series were computed to order (3^^. The number of contributing graphs, with upto 14 
bonds was 4545 for the 2-chain ladders and 5580 for the 3-chain ladders. The series are listed 
in Table IV. We use integrated first-order inhomogeneous differential approximants to 
extrapolate the series. The resulting estimates are shown in Fig. 13-14. For the susceptibil- 
ity, as a comparison, the recent Quantum Monte Carlo (QMC) results of Frischmuth et al 
p!5| and the results from our T=0 Ising expansion for 3-chain are also shown. It can be seen 
that our results agree very well with the QMC results except for the 3-chain system at very 
low temperatures. Given the recent findings that for the spin-half chain, the T = value 
is reached from finite temperatures with infinite slope |2^, one might expect the T — »• 



behavior for these 3-chain systems to be equally complex, making it very difficult to explore 
numerically. For the specific heat, our results showed good convergence upto the peak, but 
poor convergence below it. The results for 2-chain ladder are consistent with recent quantum 
transfer-matrix calculations by Troyer [O] et al. 



III. CONCLUSIONS 

We have studied the 2 and 3 chain Heisenberg-Ising ladders by a variety of different series 
expansions. Our results confirm the existence of a gap in the excitation spectrum of 2-chain 
systems, with either ferromagnetic or antiferromagnetic interchain interactions. For 3-chain 
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systems, our results arc consistent with a gapless spectrum. The T = phase diagram as 
well as the temperature dependence of the Tiniform susceptibility and the specific heat are 
also calculated. Overall, our results agree very well with previous numerical studies of these 
systems. 
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FIGURES 




0.2 0.4 0.6 0.8 1 

Xc 

FIG. 1. The phase boundary for 2-chain ladder. The asymptotic behaviour as y — — oo 
predicted by a spin-one single chain system with one-site anisotropy is also shown by the bold 
line. 
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FIG. 2. The ground-state energy per site Eq/N as function of y for 2-chain and 3-chain ladders. 
The error bars are much smaller than the symbols. 
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FIG. 3. The uniform susceptibility x at T = as function of y for the 3-chain ladder. 
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2 chain: ky=7T 
, from top to bottom 

y=2, 1.5, 1, 0.75, 0.5, 0.35, 0.1, 




FIG. 4. The dispersions of the spin-triplet excitatcd states of the 2-chain ladder with antifer- 
romagnetic interchain coupling y =2, 1.5, 1, 0.75, 0.5, 0.25, 0.1, (shown in the figure from the 
top to the bottom respectively), for ky = tt. 



11 




2 chain: ky=0 
from top to bottom 
y=0, -0.25, -0.5, -0.75, 



-1.5, -2, 
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tt/2 
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FIG. 5. The dispersions of the spin-triplet excitated states of the 2-chain ladder with ferromag- 



netic interchain coupling y =0, -0.25, -0.5, -0.75, -1, -1.5, -2, -4 , for ky 



0. 




7t/4 7t/2 37t/4 tt 

K 

FIG. 6. The dispersions of the spin-triplet excitated states of the 2-chain ladder with ferromag- 
netic interchain coupling y =-2, -1.5, -1, -0.75, -0.5, -0.25, -0.1, 0, for ky = tt. 
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FIG. 7. The dispersions of the spin-triplet excitated states of the 3-chain ladder with antifer- 
romagnetic interchain coupling y =1.5, 1, 0.75, 0.5, 0.25, 0, for ky = 0. 
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FIG. 9. The dispersions of the spin-triplet excitated states of the 3-chain ladder with anti- 
ferromagnetic interchain coupling y =2, 1.5, 1, 0.75, 0.5, 0.25, 0.1, 0, for the middle excitation 
band. 




FIG. 10. The dispersions of the spin-triplet excitated states of the 3-chain ladder with ferro- 
magnetic interchain coupling y =0, -0.1, -0.25, -0.5, -0.75, -1, -1.5, -2 for ky = 0. 
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FIG. 11. The dispersions of the spin-triplet excitated states of the 3-chain ladder with ferro- 
magnetic interchain coupling y =-2, -1.5, -1, -0.75, -0.5, -0.25, -0.1, 0, for ky = 7r/2. 




K 

FIG. 12. The dispersions of the spin-triplet excitated states of the 3-chain ladder with ferro- 
magnetic interchain coupling y =-2, -1.5, -1, -0.75, -0.5, -0.25, -0.1, 0, for the middle excitation 
band. 
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FIG. 13. Susceptibility as a function of temperature for 2-chain and 3-chain ladders from the 
high temperature series expansion, and the Ising expansion at T = (for 3-chain only), also shown 
are the QMC results of Frischmuth |15] et al as the open symbols (for 2-chain) and filled symbols 
(for 3-chain) for comparison, 
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TABLES 



TABLE I. Series coefficients for the ground-state energy per site Eq/N, the staggered/cohnear 
magnetization M, and staggered parallel susceptibility x//- Coefficients of are listed for both 
the spin-i 2-chain and 3-chain ladders with y = ±1. 



n Eq/N M X// 



2-chain ladders with y = —1 

-3/8 1/2 

2 -1.250000000000x10-1 -1.250000000000x10-1 2.500000000000x10-1 

4 -1.562500000000x10-2 -6.944444444444 xlO-^ 4.131944444444x10-1 

6 2.443214699074x10-3 1.398473668981 xlQ-^ 9.964453526878x10-2 

8 -3.391054004308x10-3 -4.293584111636x10-^ 6.561704135203x10-1 

10 -3.083929459431x10-" -9.459983497527x10-^ 2.592402875212x10-1 

12 -5.758725813964x10-" -2.517788286146X 10-^ 8.633367551243x10-1 

14 -3.508609329659x10-* -1.290819523741x10-^ 5.194241039530x10-1 

16 -3.658331257827x10-" -2.220489673071 xlO-^ 1.197591274920 

2- chain ladders with y = 1 

-3/8 1/2 

2 -1.875000000000x10-1 -1.875000000000x10-1 3.750000000000x10-1 

4 -9.114583333333x10-3 -1.176215277778x10-1 9.751157407407x10-1 

6 -6.564670138889x10-3 -1.298255449460x10-1 2.225074548290 

8 -1.011191788034x10-2 -2.567768064481x10-1 6.466470614293 

10 -4.438918908655x10-3 -3.308567310261x10-1 1.436049594755 xlOi 

12 -1.149045736404x10-2 -6.754838729347x10-1 3.661634831290 xlOi 

14 -9.014538657677x10-3 -1.059678748481 8.236161272382x101 

16 -1.860424772341x10-2 -2.087462804052 1.962459013257x102 

3- chain ladders with y = —1 

-5/12 1/2 

2 -1.111111111111x10-1 -1.018518518519x10-1 1.913580246914x10-1 

4 -9.126984126984x10-3 -2.765516082976x10-2 1.235451495119x10-1 

6 -3.430007743753x10-3 -2.229002861112x10-2 1.702591936321x10-1 

8 -4.479840217373x10-" -6.778016397592x10-3 9.228960192872x10-2 

10 -1.150888174205x10-3 -1.348552579039x10-2 1.814452046940x10-1 

12 -5.493987828590x10-" -9.549725222639x10-3 1.785680892107x10-1 

3-chain ladders with y = 1 

-5/12 1/2 

2 -1.777777777778x10-1 -1.551851851852x10-1 2.766913580247x10-1 

4 2.099353321576x10-3 -3.010635013951x10-2 2.520792079948x10-1 

6 -5.486291040711x10-3 -4.567138485811x10-2 4.781602664625x10-1 

8 -9.770232257758x10-" -2.676534614009x10-2 5.289125611753x10-1 

10 -8.540222317364x10-" -2.413110830174x10-2 6.371032432384x10-1 

12 -7.963887426932x10-" -3.026105625065x10-2 1.024475965994 
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TABLE II. Series coefficients for the perpendicular susceptibility x±- Coefficients of a;" are 
listed for 2-chain and 3-cliain ladders with y = ±1. 



n 


2-chain y = —1 




2-chain y = 1 




3-chain y = —1 




3-chain y = I 





1/3 




1/3 




11/36 




11/36 


1 


-2/9 




-1/2 




-1/6 




-13/30 


2 


1.481481481481x10" 


-2 


5.666666666667x10" 


1 


1.266534391534x10" 


-2 


4.652226631393x10 


3 


5.308641975309x10- 


- 3 


-7.143518518519x10" 


1 


-1.112134668682x10" 


-2 


-5.160474965706x10 


4 


-2.120002939447x10- 


- 2 


7.544920634921x10- 


1 


1.134331580417x10- 


- 2 


5.434191826896x10 


5 


7.283052502380x10" 


-2 


-9.203876999874x10" 


1 


-2.787971676358x10" 


-2 


-6.037870527565x10 


6 


-6.507519406167x10" 


-2 


9.542443240384x10" 


1 


3.588831815932x10" 


-2 


6.192643274088x10 


7 


-9.850515124711x10" 


-3 


-1.220670122297 




-3.496516710785x10" 


-2 


-6.538548959720x10 


8 


-3.321642365298x10" 


-3 


1.243330662897 




3.361956127668x10" 


-2 


6.836898452099x10 


9 


3.225303151489x10" 


-2 


-1.570424359869 




-4.583945317287x10" 


-2 


-7.207963780768x10 


10 


-2.164319811625x10" 


-2 


1.612928972471 




4.084885561076x10" 


-2 


7.446415246323x10 


11 


-2.229164879943x10" 


-2 


-2.181953857121 




-4.495109043570x10" 


-2 


-7.944909220697x10 


12 


-8.386607266663x10" 


-3 


2.245984866808 










13 


1.980513341011x10- 


-2 


-3.034629087270 










14 


-1.128439840082x10" 


-2 


3.158070079722 










15 


-1.625818723347x10" 


-2 


-4.546560413718 











TABLE III. Series coefficients for the dimer expansion of the 2-chain triplet spin-wave excita- 
tion spectrum e{kx,ky = tt) = y „j a„^^(l/y)"' cos(mA;2;). Nonzero coefficients a„_^ up to order 



n = 8 are listed. 



(n,m) 


0>n,m 




(n,m) 


drijiri 




(n,m) 


0>n,m 




(n,m) 


0>n,m 




(0, 0) 


1.000000000 




(5, 1) 


-2.031250000x10" 


1 


( 3, 3) 


1.250000000x10" 


1 


(5, 5) 


5.468750000x10" 


2 


(2, 0) 


7.500000000x10" 


1 


(6, 1) 


9.375000000x10" 


2 


(4, 3) 


1.250000000x10" 


1 


(6, 5) 


7.812500000x10" 


2 


(3, 0) 


3.750000000x10" 


1 


(7, 1) 


3.293457031x10" 


1 


( 5, 3) 


-9.375000000x10" 


2 


(7, 5) 


-6.042480469x10" 


2 


(4, 0) 


-2.031250000x10" 


1 


(8, 1) 


2.555847168x10" 


1 


( 6, 3) 


-3.164062500x10" 


1 


(8, 5) 


-2.657165527x10" 


1 


(5, 0) 


-6.250000000x10" 


1 


( 2, 2) 


-2.500000000x10" 


1 


( 7, 3) 


-2.222900391x10" 


1 


(6, 6) 


-4.101562500x10" 


-2 


(6, 0) 


-5.000000000x10" 


1 


( 3, 2) 


-2.500000000x10- 


1 


( 8, 3) 


2.752685547x10" 


1 


(7, 6) 


-6.835937500x10" 


2 


(7, 0) 


2.966308594x10" 


1 


( 4, 2) 


-3.125000000x10- 


2 


( 4, 4) 


-7.812500000x10" 


2 


(8, 6) 


4.957580566x10" 


2 


(8, 0) 


1.120300293 




( 5, 2) 


2.031250000x10- 


1 


( 5, 4) 


-9.375000000x10- 


2 


(7, 7) 


3.222656250x10" 


2 


(1,1) 


1.000000000 




( 6, 2) 


1.718750000x10- 


1 


( 6, 4) 


7.128906250x10- 


2 


(8, 7) 


6.152343750x10" 


2 


(3, 1) 


-2.500000000x10" 


1 


( 7, 2) 


-1.728515625x10- 


1 


( 7, 4) 


2.690429688x10- 


1 


(8, 8) 


-2.618408203x10" 


2 


(4, 1) 


-3.125000000x10" 


1 


( 8, 2) 


-5.047454834x10" 


1 


(8, 4) 


1.690521240x10" 


1 









TABLE IV. Series coefficients for high temperature series expansion of the uniform suscepti- 
bility x{T) = l3Y.iCil3'/{niT+H\), and the specific heat C{T) = J2iCiPy{ni2'+^i\). Coefficients 
Ci are listed for 2-chain and 3-chain ladders with y = 1. 



i 


X(T) for 2-chain 


X(T) for 3-chain 


C{T) for 2-chain 


C{T) for 3-chain 





8 


12 


36 


60 


1 


-12 


-20 


72 


120 


2 


12 


28 


-270 


-522 


3 


6 


-20 


-2640 


-5040 


4 


-20 


4 


90 


3270 


5 


-162 


-160 


141876 


318780 


6 


-630 


-1052 


580797 


1075767 


7 


9991 


17298 


-10663200 


-28792032 


8 


88228 


80468 


-118074186 


-291518730 


9 


-779322 


-1467200 


946669020 


3061122900 


10 


-13957358 


-12792822 


26078160405 


76820424879 


11 


55717397 


165603440 


-42521155560 


-195632449272 


12 


2827957594 


2955180058 


-6789937647207 


-22502126499801 


13 


4867299659 


-24526691326 






14 


-687967034169 


-924449102836 
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